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Abstract 



The possibility that BL Lac S5 0716+714 exhibits a linear root mean 

square (rms)-flux relation in its IntraDay Variability (IDV) is analysed. 

(-H . The results may be used as an argument in the existing debate regarding 

O , ' the source of optical IDV in Active Galactic Nuclei. 63 time series in dif- 

I , ferent optical bands were used. A linear rms-flux relation at a confidence 

y ■ level higher than 65% was recovered for less than 8% of the cases. We 

were able to check if the magnitude is log-normally distributed for eight 
!!* , timeseries and found, with a confldence > 95%, that this is not the case. 

keywords: BL Lacertae objects; individual(S5 0716-1-714) - X-Rays: 

rms-flux relation - optical: IDV. 



1 Introduction 

The linear rms-flux relation for variability in X-Rays has been established ob- 
servationally to be almost ubiquitous in a wide range of objects, starting from 
X-Ray Binaries (XRBs) to super massive black holes in the active regime (Ac- 
tive Galactic Nuclei) [37|, |36|, |32|. One of its major importance lies in the fact 



that this relation is observed even while the Power Spectral Distribution (PSD) 
shape of the light curve is stationary, implying that this relation is a funda- 
mental property of the variability process [37| . The physics behind this relation 
p\^ . is that the sources become more variable as they get brighter and it has been 

shown that it predicts a log-normal distribution of fluxes, implying that the 
variability process is non- linear [37| . |l| 

The optical/UV and X-Ray continua are thought to originate in two physi- 
cally distinct but nevertheless related regions [39j, the ogtical- X-Ray emissions 



being connected throtigh reprocessing in the disk [15|, |33| , especially in a two 



phase thermal model [34]. This framework is based on observational data and 



^Interestingly, the same linear rms-flux relation is exhibited by solar flares, even though 
the Sun does not have an accretion disk J40l |. However, the flux distribution of solar flares is 
well fltted by a power law, a feature explained by (additive) SOC (see e.g. |24| ) . while usually 
(multiplicative) processes which give a linear rms-flux also give a log-normal flux distribution. 



theoretical arguments, e.g.. the ratio of rest frame UV to X-Ray emission in- 
creases with accretion rate [17|, |39| . This framework can be relaxed and a better 
fit of the data can be obtained with multicomponent models, based on the 
working assumption that only a fraction of the optical variability is a results 
of reprocessing in the disk, while the other fraction is consistent with emission 
from a jecl (e.g. 14 1 or the energy reservoir model of |23|V 



Short time-scale variability in X-Ray may be explained by the promising 
disk-model of [21[ , which also naturally explains both the PSD of the light curve 
and the linear rms-flux relation [3|, |32| |. In fact, existence of a linear rms-flux 
relation is considered as an additional evidence for the fluctuations accretion 
scenarios because it suggests that the variability originates in the accretion 
flow 0. 

Almost ubiquitous in BL Lac objects [38| and present in other AGN flavours, 
IntraDay Variability (IDV, variability on timescales less that one day) in the 
optical regime is still an issue of debate B One of the debated aspects is the 
location of the source, where one has to take into account that microvariability 
is known to be stronger in radio loud quasars [16| . This enhanced level suggests 
that microvariability is due to processes in a jet. It was shown, for a sample of 
117 radio quiet sources, that the probability of microvariability is a smoothly 
varying function of the radio loudness ^] . Arguments of this type have led to 
simulation efforts to asses detection probabilities. Results for 10 radio quiet 
quasars, with confirmed intranight variability and with available X Ray data 
to constrain some of the features of the model, favour the assumption that the 
variability was caused by the presence of a weak jet component [7|. However, 
recent simulation results of [16| suggest that if one interprets the optical flux 
fluctuations as resulting from thermal fluctuations that are driven by an under- 
lying disk-based stochastic process, one can obtain microvariability consistent 
with observations of radio quiet quasars. 

Although there are isolated objects for which close correlation between the 
UV and X-Ray short time variability has been reported, e.g. for PKS 2155- 
304 [J, [331, or for 3C390.3 |12|, and even some low mass objects for which a 
linear rms-flux relation exists in the optical (accretin g w hite dwarfs show linear 



rms-flux in optical variability on a 10 day time scale [3^ and [ll| found optical 



linear rms flux for 3 XRBs) there is no consistent correlation between variability 



in optical-UV and X-Rray on short time scales for AGN [18| . [22| found a linear 
rms-flux relation in the U band variability on scales of tens of days, for NGC 
4151. 

We propose that a valuable argument in this debate might be offered by 



^It is interesting to note that if the response of the optical reprocessing region to incoming 
X-Ray flux is non-linear, i.e. fopt ~ fx, a large portion of the statistical properties of rapid 
(less than 100 s for object XTE J1118-I-480) variability can be explained M . 

^For variability on longer timescales, of the order of the thermal timescale, the source of 
optical/UV quasar variability may be consistently placed in the accretion disk. This conclusion 
is based on reverberation mapping arguments and on the strong correlation between both the 
characteristic timescale and rms of variability with the black hole mass (see e.g. II6II and 
references therein). 



analysing whether or not optical IDV in AGNs shows a linear rms-flux relatioqj. 

To our knowledge, there was no such systematic study so far of optical data 
from AGN that have already shown a linear rms-flux in X-Ray variability. In 
this work a first attempt is initiated, specifically for the BL Lac S5 0716+714. 
Ideally, such a study should be based on analysing simultaneous optical and X- 
Ray data of the object. This is not always possible so in the present study most 
of the optical data analysed does not have simultaneous X-Ray data counterpart. 

BL Lac S5 0716-1-714, historically identified as a radio source [301, has re- 
ceived a lot of attention in the literature and for our purposes we mention the 
work of ^, who study data series in the V, R, I optical bands, collected during 
Jan 1999 - April 2001. They analysed whether or not there is a correlation 
between the source magnitude and the amplitude of the IDV, and they found 
no such correlation. 

After a short description of the algorithm used for the analysis, we present 
the results for a number of 63 optical IDV data series for BL Lac S5 0716-1-714 
(Section [5]). Conclusions and cautionary notes follow in Section |31 

Although we focus on a limited number of datasets probing only IDV, the 
same algorithm applied to observations of longer time scales of variability might 
show very different results. There are at least two reasons for this. The first 
might be simply that IDV and longer time scale variations are of different nature. 
The second might be that the quality of long term observations is better and a 
hypothetical linear rms-flux relation in the optical could unveil. 

2 Analysis of data 

The root mean square deviation of a measured variable X is calculated as [37| 



\ 



N 



i=l 



^^^{x^-xf, (1) 



where the segment of observations is comprised of N data points, Xi is the value 
of A" at a time labelled by i and x is the medium value of X on this observational 
segment. Usually X-Ray data points are sampled with a very large frequency 
and there is a discussion whether for these data it would be better to use this 
method not on the light curve, but on its PSD [36]. For the optical data, we 
will perform this analysis directly on magnitude measurements. 

The function investigated is a; (cr„ns ) , where x is the magnitude for the optical 
data or the counts/second for the X-Ray data. The number of data points to 
mediate over so as to obtain x is 10 (a 10 data point /bin binning procedure). 



^For the timseries discussed in this work, evaluating the rms-magnitude relation serves the 
same purpose as the rms-flux. This can be seen very easily if one plots flux vs. magnitude 
for the typical parameters used in this paper. This is why throughout the body of the paper 
we consider that the magnitude-flux transformation is linear and discuss our results as if the 
behaviour of the magnitude is equivalent to that of the flux. 



2.1 Simultaneous optical and X-Ray observations 

The cases for which optical and X-Ray data were available for the same event 
are investigated first. 

Results are presented in Table [TJ noting a reference of at least one article 
where the time series has been previously analysed and where important de- 
tails such as observational setup, error analysis and weather conditions may 
be found. The Table contains information about the date when the timeseries 
were recorded (first column), the filter (second column), the slope (third col- 
umn) and intercept (fourth column) of the best linear fit and the value of the 
R^ goodness-of-fit statistics of the fit to the data (fifth column) . 



Table 1: Fits for the rms- magnitude relation for optical variability, in the case 
when simultaneous X-Ray variability data exists, and a X-Ray linear rms-flux 
relation has been found. For details see e.g. [28| for the first set, [9| for the 
second set, and [13| for the third set of data. 



Date 


Band 


slope 


intercept 


R' 




28-03-1996 
keV 


0.1-2.4 


1.169 


0.047 


0.776 


28-03-1996 


R 


-0.317 


1.004 


0.02 




5-04-2004, 
Fig.nkeV 


1.5-3.0 


1.132 


0.707 


0.915 


3-04-2004 


V 


-2.269 


1.314 


0.24 


5-04-2004 , 
Fig. [2] 


V 


0.667 


0.941 


0.491 


6 April 2004 


V 


4.752 


0.822 


0.534 




7-12-1994 
keV 


1.5-3.0 


1.668 


0.01 


0.6685 


7-12-1994 


Re 


3.873 


14 


0.215 


7-12-1994 


F98 


4.924 


13.74 


0.056 


7-12-1994 


V 


-0.733 


14.26 


0.178 



As an example for the fits, the rms-fiux spread for 5 April 2004 both in X- 
Rays and optical are shown in Figs. [1] and [2] It is clear that while for the X-Ray 
case a linear fit agrees very well with the points, for the optical case a higher 
order polynomial is needed, in this case a second order one. A cautionary note 
is necessary at this point: there are timeseries amongst those we analysed for 
which deviation from an " acceptable" linearity was due to only a finite number 
of points. This is especially obvious in Fig. [I] where if one would remove lower 
point in the plot, a linear fit would work well. 

For the date 5 April 2004, where we have almost simultaneous measurement 
for optical and X-Ray (e.g. Fig. 4 of [9|), we calculate the sample correlation 
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Figure 1: For 5 April 2004 the count rate vs. rms spread for X-Rays. 
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Figure 2: For 5 April 2004 the magnitude vs. rms spread for optical (data taken 
from [9[). In this case, the linear fit (dot-dashed) agrees very poorly with the 
data, while a second order polynomial fit (full line) produces an B? = 0.8021. 



coefficient r between the two timeseries as 



Ei=r i^^ - ^) (y^ - y) 






^T^ 



li {y^ - y) 



(2) 



where x denotes the optical timeseries, expressed in magnitude and y denotes 
the X-Ray timeseries, expressed in counts per seconco. We obtained that r = 
0.3638. A crystal clear diagnosis would have been to have \r\ G {0,1}, which 
would have meant anti-correlated (for —1), complete uncorrelated (for 0) and 
correlated samples (for l)|f|. This value of r does not provide any insight in its 



^Since we only care about variations, i.e. the behaviour of relative values, the conversion 
factors are unimportant as long as the conversion is linear. 

^Note that because magnitude varies inversely proportional to flux, the numerical value 
for the correlation coefficient for the optical flux- X-Ray flux timeseries is r = —0.3638. 



own, but in the context of our other data, it also points towards the fact that 
optical IDV does not show a linear rms-flux relation. 

2.2 Just optical data 

The rms-magnitude relation in optical IDV might prove an insightful tool for 
physical analysis of accretion in AGN. With this in mind, a summary of the 
best linear fits for a number of IDV data series is given in Tables [21 13] and ID 
The structure of the Tables is similar to that of Table [TJ 

Table 2: Fits for the rms-magnitude relation for optical variability. Details may 
be found in_J27[ for the first group, [26] for the second group, [5j] for the third 
group and 20] for the fourth group of data. 



Date 


Band 


slope 


intercept 


R' 




2-01-2000 


V 


9.709 


13.91 


0.212 


12-01-2000 


V 


-10.01 


13.54 


0.158 


25-01-2000 


V 


-2.152 


14.03 


0.016 


26 Jan 2000 


V 


3.544 


13.49 


0.266 




25-02-2003 


R 


-3.549 


13.07 


0.121 




5-03-2003 


B 


-4.213 


0.398 


0.484 


6-03-3003 


B 


0.589 


0.23 


0.074 


7-03-2003 


B 


-2.371 


0.651 


0.075 


8-03-2003 


B 


1.593 


0.043 


0.026 


9-03-2003 


B 


3.939 


-0.206 


0.154 




10-01-2007 


R 


283.6989 


-4.8039 


0.2884 


12-01-2007 


R 


-193.6 


27.12 


0.201 


23-02-2007 


R 


208.2 


2.183 


0.564 


19-03-2007 


R 


169.3 


3.354 


0.207 


20-03-2007 


R 


13.13 


13.35 


0.021 



As an example for the fits, the light curve for the B 2454826 data (Fig. [3]) 
and its rms-magnitude spread (Fig. |4]) are shown. A linear fit is superimposed, 
but it does not agree with the data. For comparison purposes, a simulated light 
curve with a log- normal fiux distribution is shown in Fig. [5] (see Section 12.31 for 
details of the simulation) . It is clear by visual inspection only that the behaviour 
of the two light curves is very different, a log-normal distribution of fiuxes leads 
to a characteristic behaviour of the peaks as compared to the mean (see Fig. (4) 



of [37[ and the discussion pertaining to that Figure) . 

A useful way to summarize the main results of the analysis performed so far 
(Tables [T]|4]) is to visualize the number of data series well fitted by various order 



Table 3: Fits for the rms- magnitude relation for optical variability, data from 



Date 


Band 


slope 


intercept 


R' 


2454824 


B 


-0.132 


1.896 


0.703 


2454826 


B 


-2.562 


14.27 


0.04 


2454828 


B 


0.01 


- 0,149 


0.7991 


2454866- 
2454867 


B 


0.088 


- 1.217 


0.253 


2454871- 
2454872 


B 


0.015 


- 0.199 


0.068 


2454872- 
2454873 


B 


-0.023 


0.359 


0.045 


2454765 


V 


-0.005 


0.081 


0.006 


2454766 


V 


-0.014 


0.201 


0.055 


2454767 


V 


-0.046 


0.633 


0.287 


2454770 


V 


0.025 


-0.348 


0.153 


2454824 


V 


-0.129 


1.787 


0.766 


2454825 


V 


0.014 


-0.194 


0.061 


2454826 


V 


-0.003 


0.05 


0.005 


2454828 


V 


0.118 


-1.657 


0.173 


2454829 


V 


0.042 


-0.579 


0.297 


2454830 


V 


-0.007 


0.101 


0.005 


2454865- 
2454866 


V 


-0.025 


0.359 


0.077 


2454866- 
2454867 


V 


0.071 


-0.952 


0.229 


2454871- 
2454872 


V 


0.043 


-0.581 


0.168 


2454872- 
2454873 


V 


-0.15 


0.225 


0.014 



degrees polynomials in a histogram plot (Fig. [§]). Only less than 8% of the data 
show a linear rms-flux relation with R^ > 65%. 

If the variability timescale r is defined as the shortest time between two 
points, xi and X2, within a time series such that \xi — X2\ > 3(7, the vari- 
ability timescales for the timeseries investigated in this paper correspond to 
T G [1.5j3.5] hours. With a value of the mass of the central object of M = 
1O*M0 [Sf and with the assumption that the perturbation was produced in 
the disk, this means that the source of perturbation is placed at a radius 
S [5.49rg, 12.81rg]. Correspondingly, we find the following values for natu- 
ral timescales of the disk |lO|: the light crossing time i/ S [1.42,3.33] hours, 
the dynamical timescale tdyn G [30.87, 110.02] hours and the thermal timescale 
itft, (ass = 1) £ [5.91,21.08] years or, more realistically, tt/i (ass = 10~^) & 



Table 4: Fits for the rms-magnitude relation for optical variability, data from 



Date 


Band 


slope 


intercept 


R' 


2454765 


R 


-0.025 


0.338 


0.12 


2454766 


R 


-0.025 


0.339 


0.197 


2454767 


R 


-0.143 


1.895 


0.82 


2454770 


R 


0.026 


0.341 


0.154 


2454824 


R 


-0.121 


1.623 


0.555 


2454825 


R 


0.064 


-0.843 


0.426 


2454826 


R 


-0.043 


0.579 


0.231 


2454828 


R 


0.074 


-1.009 


0.712 


2454829 


R 


0.021 


-0.285 


0.039 


2454830 


R 


0.042 


-0.561 


0.055 


2454865- 
2454866 


R 


-6.567 


1.38 


0.179 


2454866- 
2454867 


R 


0.084 


-1.088 


0.06 


2454871- 
2454872 


R 


0.039 


-0.513 


0.234 


2454872- 
2454873 


R 


-0.028 


0.394 


0.04 


2454824 


1 


-0.093 


0.014 


0.185 


2454825 


1 


-0.019 


0.007 


0.022 


2454826 


1 


0.017 


0.005 


0.019 


2454828 


1 


0.095 


-0.021 


0.334 


2454830 


1 


0.109 


0.0 


0.389 


2454866- 
2454867 


I 


0.041 


-0.044 


0,002 


2454871- 
2454872 


I 


0.028 


0.007 


0.14 



[591.72, 2109.03] years. Due to lack of periodicity in IDV, we were not expect- 
ing the dynamical timescale to fit the observations and it does not. The standard 
disk thermal timescale is much longer than the timescales probed in this paper 
and the standard disk viscous timescale would be even longer. However, the 
timescales of the standard disk are deduced based on a set of assumptions (see, 
e.g., [131, Chapter 7). It has been shown that relaxing the no-torque condition 
on the inner boundary of the disk causes various annuli of the disk to " communi- 
cate" on light crossing timescales [l[ and in our case the calculated light-crossing 
time fits the timescale of the observations. It is thus conceivable that the source 
of variability is placed in the disk, but not conclusive. 




0.05 0.15 0.25 

time (JD, +2454826) 
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Figure 3: the light curve of B 2454826, with error bars plotted for the magnitude 
measurements. An increment of 0.1 on the time axis corresponds to 144 minutes. 
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Figure 4: The rms-magnitude spread, with a linear fit superimposed for B 
2454826 (data taken from ^). 



2.3 Test for log-normal distribution 

A conclusive diagnosis tool for our purposes would be to test whether or not 
the optical magnitude is log-normally distributed. Statistical theory states that 
if the random variable y = log x is Gaussian distributed and if the central 
limit theorem holds, than the random variable x is log-normally distributed. A 
prerequisite for the central limit theorem to hold is that the timeseries has a 
formally infinite number of points. However, it can be shown that the conver- 
gence can be realised for a smaller number of points. Even so this condition 
is restrictive and we cannot perform this type of analysis for all the timeseries 
described in Tables [TMl 

For the purpose of this analysis, we use the chi-square goodness of fit statis- 
tics and produce Table [5] where the timeseries included have at least 100 data 
points. The first and the second column give the identification details of the 




Figure 5: Simulated light curve with the same characteristics as the B 2454826 
light curve, if the process would be log-normally distributed (see Section [2731) . 
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Figure 6: Histogram plot of the incidence of polynomial degrees which fit the 
data (according to our assumptions and to our data; the designation "7" stands 
for the cases were no polynomial fit was possible). 



time series. The third column contains the value of the x^ statistics for the 
hypothesis that the timeseries is log normally distributed, applied to the ex- 
perimental timeseries. The forth column contains the same results, but for a 
simulated timeseries. The constraint that each bin in the histogram has at least 
five members was obeyed at all times. This forced us to use only eight bins. In 
turn, for the simulations, out of the maximum eight degrees of freedom of the 
X^ distribution, three were constrained as follows: the simulated distribution 
has the same mean, variance and number of points as the observational time- 

2 
<.ref 



series. The value of the reference statistics xlef 



is 11.07. If x^ > Xr-' then 



10 



Table 5: Results for the test of the hypothesis that the timeseries are log- 
normally distributed, with Xref ~ 11-07, eight bins and five degrees of freedom. 



Date 


Band 


X"- (exp) 


X^ (sim) 


2454825 


R 


41.1184 


7.0220 


2454826 


R 


24.1224 


6.8610 


2454829 


R 


87.4848 


3.7912 


2454865 - 2454866 


R 


21.0899 


1.7466 


2454871 - 2454872 


R 


41.1184 


10.0346 


2454825 


V 


51.5006 


8.2042 


2454826 


V 


19.7716 


11.4255 


2454826 


B 


20.2235 


6.9422 



* 



histogram of the experimental data 
theoretical log-normal distribution 



13.16 13.18 13.20 13.22 13.24 13.26 13.28 13.30 



Magnitude 

Figure 7: histogram of the R 2454825 data with a log normal distribution 
superimposed. 

the tested hypothesis is rejected with a confidence of 95%. None of the time 
series we could use shows a log normal distribution. However, for the simulated 
timeseries, x^{sirn) is consistently smaller than x^T^t- For an illustration of this 
procedure, the simulated light curve of B 2454826 (Fig. [5]), the histogram of 
the R 2454825 data with a log normal distribution superimposed (Fig. [T]) and a 
simulated log normal distribution (Fig. [8]) are shown. 

3 Conclusions and cautionary notes 

The existence of a linear rms-flux relation in 63 time series in the optical data 
for the object S5 0716+714 was investigated. Although some of the data series 
show an acceptable linear fit, the general trend is that there is no such linear- rms 
flux relation in the optical band (Fig. |6|). 

We also checked if the distribution of magnitudes is log- normal (Table [S] and 
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Figure 8: simulated log normal distribution for a series with the same mean, 
variance and number of points as in Fig. [T) 



also Figs. [7]and[8|). The test was done on a restricted number of timeseries 
(eight), for reasons explained in the body of the paper. None of these timeseries 
shows a log normal distribution. 
The framework initiated by 



21 



naturally explains this relation in the X- 
Ray band by assuming independent variations in the accretion rate at different 
radii and it is, by its assumptions, a disk model. If one agrees to extend this 
argument to say that all disk variability should then exhibit a linear rms-flux 
relation, one can say that the optical IDV variability analysed in this work is 
not of disk nature. 

However, it is obvious that this preliminary result should be taken with 
the necessary cautionary notes. The rms-magnitude behaviour might prove to 
be very different for longer timescale variations. It is thus necessary that all 
groups of timescales be analysed in our proposed framework. Because of the 
binning procedure, a better statistics would be obtained if the length of each 
individual timeseries would be a generous multiple of 10. The probability of 
IDV detection itself increases with observation time (see, e.g., [31j). Longer 
timeseries are also needed for the proper testing for a log-normal distribution 
of the magnitudes/fluxes. In order to confirm or disprove our initial analysis 
for BL Lac S5 0716+714 and to add a greater degree of generality to the result 
by including other objects in the analysis, all of these notes must be taken into 
account in the future. 
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